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PREFACE 

The image analysis phase of this investigation vas designed to use 
the time-lapse analysis capabilities of ESIAC (Electronic Satellite 
Image Analysis Console) at the Stanford Research Institute. This equip- 
ment was being designed and tested under MSA Contract MS 5-21841 con- 
currently V7ith this investigation. In addition, this investigation used 
the photographic extractions of basin characteristics being produced 
from satellite imagery by the Topographic Division, U.S. Geological 
Survey, as part of a thematic mapping investigation under MSA Coiitract 
S-70243-AG. 

The objective of the investigation was to test two hypotheses. 

First, it was postulated that indices of dynamic basin characteristics 
could be extracted from imagery acquired by the Earth, Resources Technology 
Satellite (now. Landsat) . Second, it was postulated that these indices 
could be used quantitatively to improve equations for estimating stream- 
flow characteristics. 

In order to test the hypothesis that basin characteristics derived 
from. Landsat imagery can be used to improve equations for estimating stream- 
flow, it was necessary first to extract the data from the imagery and 
then to use the data in a simple analytic experiment. Several techniques 
were tried for accurately extracting and measuring basin characteristics; 
some were representative of state-of-the-art at the time they were tried. 
The analytic experiment involved computing and comparing two sets of 
multiple linear regression equations, where one set contained basin char- 
acteristics derived from Landsat imagery. This analytic experiment was con 
sidcred the primary objective limiting the scope of work. 


• It. is concluded that some hasin characteristics can be extracted 
easily' from a single set of Landsat imagery with the aid of multispectral 
image analysis systems that employ film-density discrimination techniques 
only. Two of the easily-extracted basin characteristics can be used to 
substantially impro've the accuracy of equations used to estimate 12 out 
of 40 • streamf low characteristics for streams on the Delmarve Peninsula of 
Delaware, Maryland, and Virginia. This improved accuracy could allow 
a reduction in effort needed to collect data that is used to define these 
streamflow characteristics at gaged sites. The savings in manpower and 
money could be applied to meeting other current and anticipated data 
needs . 

Mew image analysis techniques need to be tested as to their capabil- 
ities for operational output of basin characteristics data. Techniques 
using digital data from computer compatible tapes are a likely source of 
data for expanding the matrices used to develop better equations for 
estimating streamflowi 

Consideration might he given to testing basin characteristics derived 
from satellite data for improving estimates of streamfloxir in other phy- 
siographic regions in the nation. 

Consideration might also be given to the possibility of applying 
results of this investigation nationally by developing an expanded basin 
characteristics file for correlation with the national streamflow data 
files. 
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Multiply by 

Metric 

ft^/s 
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mi (square miles) 
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Icm^ (square kilometres) 
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niPROVING ESTIMATES OF STREMO-'LOW CHARACTERISTICS 
USING LANBSAT-1 (ERTS-1) IMAGERY ■ 

By Este F, Hollyday, Nashville, TN. 

Work done in cooperation with the National Aeronautics 
and Space Administration 

Abstract , — Imagery from the first Earth Resources Technology 
Satellite (renamed Landsat-1) was used to discriminate physical features 
of drainage Lasins in an effort to improve equations used to estimate 
streamflow characteristics at gaged and ungaged sites. Records of 20 
gaged basins in the Delmarva Peninsula of Maryland, Delaware, and 
Virginia were analyzed for 40 statistical streamflow characteristics. 
Equations relating these characteristics to hasin characteristics were 
obtained by a technique of multiple linear regression. A control group 
of equations contains basin characteristics derived from maps. An exper- 
imental group of equations contains basin characteristics derived from 
both maps and imagery. Characteristics from imagery were forest, riparian 
(stream-bank) vegetation, water, and combined agricultural and urban land 
use. These basin characteristics were isolated photographically by tech- 
niques of film-density discrimination. The area of each characteristic 
in each basin was measured photometrically. 
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Comparison of equations in the control group with corresponding equa- 
tions In the experimental group reveals that for 12 out of 40 equations 
the standard error of estimate was' reduced by more than 10 percent. As 
an example,, the standard error of estimate of the equation for the 5-year 
recurrence-interval flood peak was reduced from 46 to 32 percent. Sim- 
ilarly, the standard error of the equation for the mean monthly flow for 
September was reduced from 32 to 24 percent, the standard error for the 
7-day, 2-year recurrence low flow as reduced from 136 to 102 percent, and 
the standard error for the 3-day, 2-year flood volume tras reduced from 30 
to 12 percent. 

It is concluded that data from Landsat imagery can substantially 
improve the accuracy of estimates of some streamf low characteristics 
at sites in the Delmarva Peninsula. 

INTRODUCTION 

Nationally the U.S. Geological Survey operates more than 10,000 
stations that are used for gaging streamfiow. Records of stage collected 
continuously at gaging stations are converted to streamfiow and are pub- 
lished annually. These data are useful to planners, engineers, and water 
managers for designing water supply reservoirs, controlling pollution, 
designing bridges, managing flood plains, forecasting and managing floods, 
producing power, and designing and maintaining navigation facilities. 

The purpose of this investigation is to determine if imagery from Landsat 
can improve on one aspect of this program, namely regionalization of 
streamfiow information, or the transfer of streamfiow records from gaged 
to ungaged sites. 
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In 1970 the Survey initiated a study to evaluate the national atreamflow 
data program which had evolved over the previous 80 years (Benson and 
Carter, 1973). Existing data were evaluated in terms of newly-established 
goals, and a proposal for continuation of a revised program was released 
to the open file (Forrest and Walker, 1970). 


- 3 - 



As part of this evaluation, the records for each gage on unregulated 
streams were analyzed to derive. statistical measures- of flow termed 
strearaflow characteristics, which include average flows, variability in 
average flows, flood peaks, flood volumes, and low flows. Using multiple 
linear regression techniques the streamflow characteristics were then 
correlated with basin characteristics, which are selected physiographic 
and climatic features of the corresponding drainage basins. This regres- 
sion analysis was done in part to transfer the streamflow records from 
gaged to ungaged sites, and in part to evaluate the record. The statis- 
tical model that was used is: 

b(l) b(2) b(n), 

Y = aX(l) X(2) ...X(n) 

where Y is a streamflow characteristic, X(l) to X(n) are basin char- 
acteristics, and a, b(l) to b(n) are the regression constant and 
coefficients. 

The equations so generated were used to determine data collection 
needs by comparing the accuracy of the equations with the accuracy goals 
specified for estimates of flow characteristics at ungaged sites. These 
accuracy goals were given in terras of equivalent years of record. This 
means it was specified that information provided for any ungaged point 
on a stream should be equivalent in accuracy to that which would have 
been attained by an actual record of a selected number of years (10 or 25) 
at that point. Accuracy goals in terms of equivalent years of record 
in a given stS.te or region were converted to standard error in percent 
of mean using the methods described by Hardison (1969). Independent 
of specifying these accuracy goals, a value called the standard error of 
estimate, was computed for each regression equation. 
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This value is a general index o£ the accuracy of estimates obtained by 
use of the equation. By comparing the standard error of estimate of 
each streamflow characteristic with the accuracy goal for that char- 
acteristic,. it was possible to judge the degree to which some goals 
had already been achieved by existing data and to judge the need for con- 
tinued data collection. According to the concepts used in the eval- 
uation study (Benson and Carter, 1973) , about 10 percent of the present 
strearaflow program effort should he redirected to areas of higher 
priority. 



Since 1970 the streamflow' data program is, being reevaluated period- 
ically, and this investigation is pertinent td the continuing evaluation 
effort in that it tests the impact of additional basin characteristics 
on estimates of streamflow characteristics. If basin characteristics 
derived from Landsat data are added to an equation and if this addition 
results in a reduced standard error of estimate of the streamflow char- 
acteristic then the neX'f standard error may be less than the accuracy goal. 
In such a case, a reduction in effort to collect data on that character- 
istic would be in order. In this way, additional basin characteristics 
derived from satellite data could have a substantial impact on the stream- 
flow information program. 

The basin characteristics used in the regression analysis that was 
part of the evaluation were compiled from U.S. Geological Survey topograp- 
ic maps and from National Weather Service Climatological Data. The maps 
are not the most suitable source of information on basin characteristics. 
The maps Vary in scale and detail; they also vary in age of photography 
used for compiling land-cover inforraatioB.. In order to have maximum 
ground visibility, mapping photography is usually taken in early spring 
before lea£ing-out of trees which generally coincides with conditions of 
high water in the eastern half of the nation. The scale is commonly 
1:20,000 or about 3 mi (5 km) on a photograph edge. 
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On July 23j 1972, the D.S. National Aeronautics and Space Admin- 
istration (NASA) launched the first Earth Resources' Technology Satellite 
(renamed Landsat-1) capable of repeatedly and uniformly imaging the 
Earth. The opportunity is presented. to evaluate 'basin characteristics 
extracted from satellite images. Although the resolution of the Landsat-1 
system cannot compare with standard mapping photography, several advan- 
tages are foreseen in using Landsat imagery as a source of basin char- 
acteristics rather than topographic maps or the photographs from x^hich 
they were compiled. Landsat can provide seasonal information on land- 
cover conditions rather than early spring conditions only. A single 
Landsat image covers an area of 13,200 rai^ (34,200 km^) or 115 mi (185 km) 
on an image edge. The viex^ has nearly uniform lighting conditions allowing 
more uniform application of criteria for extracting basin characteristics 
than a photographic mosaic of comparable coverage. An interpreter can 
select those characteristics believed to be most closely related to hydrol- 
ogy rather than accepting the standard cultural information on topographic 
maps. 

Many previous applications of Landsat imagery to x/ater problems have 
emphasized the inventory or mapping of land-cover types or conditions 
that should be related to hydrology. This investigation sought to test 
quantitatively the inferred relationships between streamflow and selected 
land-cover types. This report describes results of testing the hypotheses 
that hydrogically-signif leant basin characteristics can be extracted from 
Landsat-1 imagery, and that these characteristics can he used quantitative- 
ly to improve equations for estimating streamflow characteristics. 
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PHYSIOGRAPHY OF STUDY A31EA 

The study area selected for testing Landsat imagery in this investi- 
gation encompasses most of the Delmarva Peninsula and part of the adjacent 
mainland covered by Landsat-1 image 1079“ 15 133 and succeeding images with 
comparable coverage (fig , 1) . The area lies within the Central Atlantic 

Figure 1, — (caption on next page) belongs near here* 

Regional Ecological Test Site (CARETS) of the U.S. Geological Survey. 

It is part of the Chesapeake Bay group of Landsat investigations funded 
by the National Aeronautics and Space Administration. 

Twenty gaged drainage basins were selected within the study area 
(table 1). Drainage areas range from 3.85 rai^ (9.97 km^) to 113 mi^ 

(293 km^) and average 24.6 mi^ (63.7 lon^) . From 11 to 31 years of 
streamflow data were available for each gage when the data were eval- 
uated in 1970. 

Mean annual precipitation ranges from 46 in (1170 mm) to 48 in 
(1220 mm). Mean annual temperature ranges from 12 ° to 14° C, and ex- 
tremes are moderated by neai'by Chesapeake Bay and Atlantic Ocean, 

During winter there is a 45 percent probability that the area will be 
cloud covered during a satellite overpass. By using parts from as many 
as five successive images of the same scene, however, complete coverage 
each season is practically assured. 
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used In multiple regression analysis. 


Drainage 

area» Period of 

•in m*r2 ,-jann'r^ nnnli 


Longitude 

in 

record analyzed 

0754010 

3.85 

1956-67; Annual max. 1952-56 

0753757 

9.35. 

1933; 1943-57; Annual max. 1958-67 

0751939 

7. OB 

1956-67 

0732031 . 

5.24 

1943-67 

0751930 

60.5 

1949-67 

0752819 

44.9 

1349-67 

0754018 

5.8 , 

1951-67 

0753411 

19.5 

1936-67 

0753341 

75.4 

1943-67 

.0752858 

16.7 ■ 

1951-67 

0754024 

44.8 

1943-67 

0754734 

7.10 

1950- 

0755250 

15.0 

1951- 

0754709 

113 

1948- 

0755815 

5.85 

1950- 

0760630 

8. 09 

1951-56; Annual max. 1957-67 

0755140 

22.3 

1948-67 

0760054 

10.5 

1951-67 

0755851 

12.5 

1951-56; Annual max, 1957-65 

0760610 

5.31 

1948-58; Annual max, 1965-67 




The basin for station number 01579000 lies within the Piedmont 
Plateau physiographic province; all other basins are situated on the 
Peninsula within the Atlantic Coastal Plain province^ This part of the 
Coastal Plain is a flat, low almost featureless plain. Maximum eleva- 
tions throughout the Peninsula rarely exceed 80 feet (24 m) above mean 
sea level. Maximum elevations in the Piedmont may exceed 400 feet (120 ra) 
near the north end of the Peninsula. Relief within any one square mile 
may have a value as much as half the local maximum elevation. The shore- 
line of the Peninsula is extremely broken and sinuous along the 
Chesapeake Bay and is characterized by very small relief and tidal 
marshes. The shore of the Atlantic Ocean is composed of a long line of 
barrier beaches with lagoons on the landward side. All gaging stations 
are located away from the shore and above tidal influence. The Coastal 
Plain is underlain by a series of southeasterly dipping layers of uncon- 
solidated sand and clay with a subordinate amount , of gr'avel. These 
layers overlie the weathered crystalline rock of the Piedmont and thicken 
to the southeast from the northwest edge of the Coastal Plain. These 
materials are drained hy sluggish rivers, many of which have been or 
are being channelized above tide water for the purpose of improving 
the drainage of agricultural land. 



Before colonial time, the area was completely covered with forest. 
The amount of land cleared for agriculture has varied throughout history. 
Currently, second- and third-growth forest is largely restricted' to 
river flood plains, swamps, and wet upland depressions in the northern 
two-thirds of the area and also to state forest preserves in the southern 
one- third of the area. Upland forest is composed predominantly of pine, 
mostly Virginia pine in the northern part of the Peninsula and loblolly 
pine in the south. Lowland areas have a distinctive swamp hardwood 
forest composed of pin, willow, and swamp oak, red and black gum, red 
maple, river birch, yellow poplar, sycamore, beech, and walnut (Yokes, 
1957). Permanently-flooded or tidal, fresh-water areas in the south 
have extensive stands of bald cypress. Agricultural land in the northern 
half of the study area, is used predominantly for dairy farming, live- 
stock, and feed grain, particularly corn. . In the southern half it is 
used for truck farming and poultry production. 


Level 1 land use ([Anderson and others, 1972) within the 20 selected 
drainage basins was predoniinantl)f forest and agriculture in 1970 with 
the latter category being most common in the northern half of the study 
area (E, J. Pluhowslci, written commun. 1974). 


Level 1 land-use 

Land use within basins. 

In percent 

category 

range 

average 

Urban and built up 

0.0 - 5.1 

0,5 

Agriculture. 

20.2 - 96.2 

58.4 

Eorest 

3.8 - 79.6 

40,9 

Water 

0.0 - 0.8 

0.1 


LANDSAT-1 

Landsat-1 is an experimental satellite for demonstrating that remote 
sensing from space with unmanned Satellites is feasible and can provide 
valuable data to assist efficient management of water and other earth 
resources (National Aeronautics and Space Administration, 1972). 

Although the satellite was launched in July 1972 with a design life of 
1 year, it has continued to provide imagery and to relay data through 
June 1975. In order, to allow systematic, repetitive imaging of the 
Earth under nearly constant lighting conditions, it was launched and 
maneuvered into a circular, near-polar orbit that allows the satellite to 
repeat its ground trace at the same local solar time every IS days. 

This is 0942 hours at the equator and about 10 minutes earlier for most 
of the lower 48 states and Hawaii. 



The. satellite contains two ■ multispectral Imaging systems j a data 
relay system, and support systems. The Multispectral Scanner (MSS) is 
a line scanning device with arrayed detectors rather than film or sensi- 
tiaed-phosphor plates. It detects daylight solar energy reflected from 
the Earth's surface in the visible .and-near (non-thermal) infrared region 
of the spectrum. This energy passes through a single optical system 
that allows it to be recorded in four spatially-registered spectral bands. 


Multispectral Scanner 
hand 

4 

5 

6 
7 


Wavelength, in 
micrometres 
0.5 - 0.6 
0.6 - 0.7 
0.7 - 0.8 
0.8 - 1.1 


The human eye is sensitive to an 0.4 to 0.7 micrometre wavelength hand. 
The satellite system scans from west to east at right- angles to the path 
of travel. Forward motion of the satellite provides continuous coverage 
along the orbital track. The detector response is sampled, bit encoded. 


and transmitted to a ground station. During image processing, the con- 
tinuous data stream is located geographically, corrected geometrically 
and radiometrically, and framed to produce a 55-mm image on 70-mni film. 
Each frame covers an area of 34,200 Ion . Further processing provides a 
variety of products and formats including black and white or color-com- 
posite images with a scale of 1:1,000,000 on 9-l/2~inch C240-inm) film. 
Each image has 10-percent forward lap with immediately previous and 

succeeding images and about 20 percent side-lap for imagery of most of 
the nation. Spatial resolution for the images averages about 80 m, but 
may be considerably higher or lower depending upon scene contrast.; 
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BASIN . CHARACTERISTICS 
CHARACTERISTICS FROM LANDSAT IMAGERY 
Spectral Reflectance 

A basin characteristic observed in Landsat itnagcry absorbs and re- 
flects solar radiation in its o\m. way. The incoming solar radiation 
striding an object may be compared with the resulting outgoing or re- 
flected radiation for small increments of wavelength. The resulting 
data can be used to produce spectral reflectance curves for the selected 
object (fig. 2). Dry snow reflects greater than 50 percent and clear 


Figure 2. — (caption on next page) belongs near here. 

water reflects less than 10 percent of solar radiation over the band 
xfidth 0.3 to 1.1 micrometres. As a result dry snow appears very bright 
in cloud-free imagery for all four MSS bands, and water appears very 
dark, especially in bands 6 and 7. Color films that have emulsion 
layers sensitive to blue, green, and red light in the visible spectrum 
portray forest as green because maximum reflected solar radiation occurs 
in the green part of the visible spectrum. An infrared color composite 
of Landsat bands 4, 5, and 7, however, portrays forest as red because 
maximum reflected radiation occurs in the near-infrared or non-visible 
part of the spectrum beti^reen 0.8 and 1.1 micrometres, band 7, 






Figure 2. — ^Generalized spectral reflectance curves for five land- 
cover categories compared to Landsat spectral bands. 
Curves for concrete, asphalt, forest, and water are 
modified from Root and Miller (1971); for dry snow 
from Serebreny and others (1974). 
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In theory, a specific basin characteristic can be isolated in 
Landsat imagery from others by .comparing film densities (which are re- 
lated to spectral reflectance); however, difficulties are frequently 
encountered when attempts are made to isolate a characteristic using 
film densities only. For a variety of reasons there may be significant 
overlap among spectral curves. For example, a thin cloud cover between 
ground and satellite tends to increase the apparent spectral reflectance 
of all terrain features. In this way water with a thin cloud cover may 
easily be confused with asphalt in another part of the same .scene where 
there is no cloud cover. Suspended sediment also increases the spectral 
reflectance of water. Thus water with sediment may be confused with 
asphalt in coastal cities even in cloud-free imagery. Conversely, shadow 
from either cloud or terrain tends to decrease the apparent spectral reflec 
tance. In this way, asphalt in shadow may be confused with water in open 
sunlight. In addition to these problems of isolating a specific basin 
characteristic in a single Landsat image, seasonally, the angle of solar 
illumination changes, trees lose or gain their leaves, and soil moisture 
varies. Accordingly, criteria for isolating basin characteristics must 
also change from scene to scene throughout the year. 

These difficulties limit the success of isolating a basin character- 
istic that has been selected solely on the basis of hydrologic signifi- 
cance. Some preselected characteristics do not have a unique spectral 
response that is always detectable with the Landsat system. As a result 
the selection of basin characteristics to be derived from Landsat imagery 
is a process involving compromise between hydrologic significance and 
ease of extraction. 
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selected Basia Characteristics 

Four basin characteristics were finally selected after giving con- 
sideration to hydrologic significance, ease of extraction,., and availabil- 
ity of hydrologic and remote-sensing data in the study area* The four 
characteristics and their symbol designations are:, forest (Uf) , water 
(Uw ) 3 riparian vegetation (Urv) , and combined agricultural arid urban land 
use (Uau) ... 

Forest (Uf) is the relative area of a drainage basin that is covered 
by trees, expressed in percentage of total basin area. In the study area 
this includes upland pine forest as well as lowland, swamp hardwood forest. 
Forest appears to affect streamflow indirectly by affecting rates of 
evapotranspiration, precipitation interception, and the accumulation and 
melting of snow. Forest cover was found significant at no less than the 
5 percent level in estimating 19 out of 41 streamflow characteristics 
(Forrest and Walker, 1970). Forest has a unique spectral response detect- 
able with the Landsat system. 

Water (Uw) is the relative area of a drainage basin that is covered 
by a water surface visible from iSibove, expressed in percentage of total 
basin area. In the study area it is the area of mill ponds, industrial 
waste lagoons, and flooded fields. Water is a measure of surface water 
storage which affects evaporation, peak flows for short- and intermediate- 
recurrence intervals, and low flows of natural streams. Storage was 
found significant at no less than the 5 percent level in estimating 2 
out-. of 41 streamflow characteristics (Forrest and Walker, 1970). Clear 
water has a unique spectral response that is easily detected with the 
Landsat system. ; ' 
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Riparian vegetation (Urv)' is the relative area of a drainage basin 
that Is covered by vegetation situated on or near the bank of a stream 
or other water body, expressed iii percentage of total basin area. 

In the study area it is area of shrubs and trees in swamps, marshes, 
and seasonally-flooded depressions. The roots of riparian vegetation 
have constant access to water and thereby affect streamflow, particularly 
intermediate and low flow, by affecting rates of evapo transpiration. 

The roots can either withdraw water , directly from a stream or intercept 
water moving through the ground to a stream. Riparian vegetation was 
not used previously in estimating streamflow characteristics. In some 
areas riparian vegetation has a unique spectral response, that can be 
detected in Landsat imagery for winter months. 

Combined agricultural and urban land use (Uau) is fields, pasture, 
buildings, roads, and sand pits, expressed in percentage of total basin 
area. In the study area it is all areas not covered by either forest 
or water. Wo attempt was made to extract this basin characteristic 
from Landsat imagery. It was arbitrarily measured as the difference be- 
tween total basin area and the combined area of forest (Uf) and water 
(Uw)i 
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Image Analysis 

To make a quantitative test of the inferred relationships bett^een 
strearaflow and selected basin characteristics derived from Landsat 
imagery the basin characteristics are first measured and then used in an 
analytxc experiment. Figure 3 shows the key steps in this analysis. , 


Figure 3.- — -(caption on next page) belongs near here. 

Cloud-free imagery for each season is selected from catalogued Landsat 
data. The imagery is inspected to estimate the ease of isolating the 
selected characteristics from a single band, or a combination of bands, 
over one or more seasons. Characteristics such as forest are then 
isolated and extracted. The relative area of each characteristic in 
each basin is then measured. After Image analysis, the characteristics 
derived from Landsat data are merged with an available matrix of map- 
derived basin characteristics. The new matrix is then correlated with 
a matrix of streamflow characteristics by a technique of multiple re- 
gression. 

During image analysis, effort was concentrated within the 20 gaged 
basins. A mask was used to black-out all of an image except the areas 
within the basins. The basin mask was prepared by transferring the 
basin boundaries from topographic maps to opaque film. The opaque layer 
of this film was then removed from the area enclosed by the boundaries. 
Geographic control points were then added to allow registration of the 
mask with imagery having a scale of 1:1,000,000. 
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APPLICATION IN 
ESTIMATING STREAMFLOW 


Flow diagram of technique for improving equations used 
to estimate streamflow characteristics. Forest areas 
are extracted from Landsat image 1079-15133, are meas 
ured basin by basin, and are added to characteristics 
used previously in developing regression equations. 


Reproducibility of thf 
rtftmiKAT. PAGF IS POOR 


tianual extractions of Eotest and water were prepared as guides 
to automated image analysis using multispectral image analysis systems 
and using photography. A positive film transparency of band 5 of image 
1079-15133 was overlain with the basin mask and a clear film was placed 
over the mask. The distribution of the selected basin characteristic 
was estimated visually and dra'im on the clear film by subjectively 
applying criteria of film density, texture, shape, and terrain posi- 
tion relative to the drainage pattern. Areas believed to meet the 
criteria were inked- in to complete the extraction. 

Seasonal changes that might be useful in identifying and extract- 
ing basin characteristics were detected by time-lapse processing of two 
or more images of the same scene (Serebreny and others, 1974). Color 
displays made of a magnified part of image 1079-15133 of October 10, 

1972, and of image 1187-15140 of January 26, 1973, using ESIAC 
(Electronic Satellite Image Analysis Console at Stanford Research 
Institute) revealed that riparian vegetation could be distinguished from 
upland hardwood forest in the winter image of basin 01485000. Values 
of scene radiance along a cursor trace through the image of the basin 
indicated that the variation in film density (directly related to scene 
radiance) was just large enough to isolate riparian vegetation by equi- 
densitometry . ' 


. Photomechanical extractions were made of. forest, riparian veg- 
etation, and xsrater using one or more bands of images 1079-15133 and 
1205-15141 of February 13, 1973 1 .' The film density of a selected basin 
characteristic throughout the image was measured along with the density 
of selected steps in the image gray scale. A suitable range in density 
X'fas selected for isolating the characteristic and was specified for 
photochemical control in a' high-quality graphic arts photographic lab- 
oratory. The laboratory used suitable density- isolating films to pro- 
duce binary extractions in the form of photographic transparencies. 

Two or more of these transparencies can form a composite or sandwich 
that reduces unwanted or spurious data while isolating the desired 
basin characteristic. The photomechanical extractions used in this 
investigation were produced in the laboratory that is part of the 
Autographic Theme Extraction System of the Earth Resources Observation 

Systems program of the Department of Interior (Smith, 1973) . 

Measurement of the relative area of each characteristic in a gaged 
basin was made photometrically. Light from a photographic enlarger was 
focused upon an integrating photometer. The basin mask was introduced 
between the light source and photometer. The total illumination pass- 
ing through the transparent area inside a basin boundary on the mask is 
a measure of the total area of the basin. A photomechanical extraction 
with transparent areas representing the basin characteristic was then 
registered with the mask. The total illumination passing through both 
basin mask and extraction is a measure of the area of the basin character 
istlc only. The ratio of illumination measurements was expressed as 
percent of basin covered by the selected basin characteristic. The 
photometric equipment was calibrated by introducing circular holes of 
known diameter between the light source and the photometer. 
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The ratio of xll-uminations for different hole diameters agreed within 
5 percent with the ratio of areas calculated for the corresponding 
holes. The equipment selection and setup were, conceived by, W. E. Evans, 
Stanford Research Institute. Basin' characteristics determined by this 
method are summarized in table 2. 







Table 2. — Basin characteristics derived from Landsat-1 imagery. 


Area of basin covered by given characteristic, 
in percent of total basin area 


Station 

Forest 

Water . 

Riparian 

Combined 

number 

(Of) 


vegetation 

agr, and urban 

CUw) 

(Urv) 

(Uau) 

01483200 

43.6 

1.1 

52.8 

55.3 

01483500 

14.2 

0.0 

21.0 

85.8 

01484300 

29.6 

1.8 

6.1 

68.6 

01484500 

23.9 

0.0 

3.9 

76,1 

01485000 

40.8 

0.0 

11.9 

59.2 

01485500 

79.4 

0.0 

5.4 

20.6 

01486000 

67.4 

0.0 

10.8 

32.6 

01486500 

39*4 

0.7 

0.8 

59.9 

01487000 

29.3 

0,0 

3.5 

70.7 

01487500 

■ 59.7 

1.0 

2.4 

39.3 

01488500 

32.3 

0.0 

18.0 

67.7 

01489000 

19.5 

0.0 

3.0 

80,5 

01490000 

55.2 

0.9 

3.6 

43.9 

0 1491000 

40.4 

. 0.4 

39.1 

59.2 

01492000 

36.2 

0.0 

3.0 

63.8 

01492500 

31.3 

0.0 

13.8 

68.7 

01493000 

23.3 

0.8 

35.6 

75.9 

01493500 

6,0 

0.0 

6.7 

94.0 

01494000 

20.2 

0.0 

18.0 

79.8 

01579000 

21.3 

0.0 

32.2 

78.7 



CHARACTERISTICS FROM MAPS AND CLIMATOLOGICAL RECORDS 

The twelve physiographic and climatic characteristics used in the 
regression analysis that was part of the evaluation for Maryland and 
Delaware (Forrest and Walker, 1970) include the following: 

A, drainage area, in square miles, contributing to surface runoff, 
derived from topographic maps and shown in the latest U.S. 
Geological Survey streamflow report, 

S, main channel slope, in feet per mile computed by the 85- to 
10-percent method (Benson, 1962), 

L, stream length, in miles, measured along the main channel from 
gage to basin divide, . 

E, mean basin elevation , in feet above mean sea level, measured 
from topographic maps by the grid sampling method (20 to 80 
points in basin were sampled) , 

St, area of lakes, ponds, and swamps, in percent of contributing 
drainage area, measured by the grid sampling method, 

F, forest area, in percent of contributing drainage area, meas- 
ured by the grid sampling methods, 

Si, soil index, a relative measure of potential maximum infiltra- 
tion capacity in- Inches, estimated from data provided by the 
U.S. Soil Conservation Service, 

P, mean annual precipitation in inches, from The National Weather 
Service, "Climates of States;" grid sampling method used on 
isohyetal maps. 



124 j 2, precipitation intensity; TnaximuiK 24-hour rainfall, in 
inches, expected on the average of once every 2 years; 
estiioated from U.S. Weather Bureau Technical Paper 29, 

Sn, mean annual snowfall, in inches, from The Rational Weather 
Service, "Climates of States," 

Tl, mean minimum' January temperature, in degrees Fahrenheit, 

from The National Weather Service, "Climates of States," and 
T7, mean maximum July temperature, in degrees Fahrenheit, from 
The National Weather Service, "Climates of States." 

Values for these twelve basin characteristics for each of the 
20 gaged basins used in this investigation are available from the 
U.S. Geological Survey computer file, Streamf low/Basin Character- 
istics, and are published in previous reports (Forrest and Walker, 
1970) . The hydrologic significance of each basin characteristic 
is discussed in previous reports (Benson, 1962; Thomas and Benson, 
1970). Identical values were used in this investigation and in 
the original streamf low data evaluation (Forrest and Walker, 1970). 



REGRESSION ANALYSIS 


Relations between streamflow characteristics (dependent variables) 
and drainage basin characteristics (independent variables) were derived 
by multiple linear regression analysis. This analysis, provides an equa- 
tion of the statistical relation between a streamflow characteristic, 
and selected basin characteristics. It also provides a measure of the 
accuracy of the relation defined for the sample population (knoim as 
the standard error of the estimate of the dependent variable, or simply 
standard error of estimate) . The standard error of estimate is a meas- 
ure of the spread of the data about the line of relation. It is a sta- 
tistical parameter such that the value of the streamflow characteristic 
estimated mth the equation for 2 out of 3 gaging stations will, on the 
average, plot within one standard error of the curve of relation. In 
addition the value for 19 out of 20 gaging stations will plot within 
two standard errors of the curve. 


- 28 - 


equations are calculated for each of the 40 streamflow character- 
istics using the 20 gaged basins in the study area and 12 basin character- 
istics derived from maps only, then these equations are comparable to the 
equations derived in streamflow program evaluation for Maryland and 
Delaware (Forrest and Walker, 1970) and constitute a control group for 
experimentation. The standard error of estimate will provide a measure 
of the accuracy of the equation. If equations are then calculated for 
each of the 40 streamflow characteristics from the same 20 gaged basins, 
the same 12 map-derived basin characteristics, and 4 basin character- 
istics derived from Landsat imagery in the same analysis, then these 
equations constitute an experimental group for testing characteristics 
derived from Landsat. The new standard error of estimate will provide 
a measure of the accuracy of the new equation. Any reduction in the 
standard error of estimate between a new or experimental equation and 
an old ox control equation represents an increase in the accuracy of the 
estimate of that particular streamflow characteristic. The improvement 
is due solely to including basin characteristics derived from Landsat 
imagery. 


- 29 - 


Calculations required for the stepforward regression analysis were 
performed by digital computer using a group of computer programs known 
as STATPAC. The program group eliminated any indefinite values from the 
dependent variables, added a very small constant (0.0001) to those var- 
iables which might be expected to go to zero and transformed all dependent 
variables and four independent variables to their logarithms. Program 
D0094, Multiple Linear Regression (stepforward) , first provided the follow- 
ing statistical parameters', means, standard deviation, and correlation 
matrices of all variables. 

The Stepforward program computed a set of equations by starting with 
the most effective independent variable then adding the next most effec- 
tive variable and additional variables until the accuracy of the equation 
was not significantly improved by any additional variables. After com- 
puting each equation it provided several parameters including regression 
constant and coefficients, multiple correlation coefficient, standard 
error of estimate of the dependent variable, and percent of the total sums 
of squares of the dependent variable that are explained. In addition, 
the program tabulated the observed, the calculated, and the residual of 
each streamflow characteristic (dependent variable) for each of the 20 
basins. 


The streamflow characteristics (dependent variable) used in the 
regression analysis for the streamflow data evaluation for Maryland and 
Delaware (Forrest and Walker, 1970) were defined at, 105 gaging stations 
and include the full range of flow. In computing these characteristics, 
frequency curves were not extrapolated beyond twice the length of 
record. These characteristics Include the following: 

FT, annual flood peak, in cubic feet per second, of T-year recur- 
rence interval; the recurrence intervals of 2, 5, 10, 25, and 
50 years. are denoted in this report as P2, P5, PIO, P25, and 
P50 respectively, 

QA, mean annual discharge, in cubic feet per second, defined as 
the mean of the annual means, 

SDA, standard deviation of mean annual discharge, in cubic feet per 
second, 

QM, mean discharge for the M-calendar-month, in cubic feet per 

second; the M refers to the numerical order of the month begin- 
ning with January as 1, 

SDM, standard deviation of mean discharge for M-calendar-month, in 
cubic feet per second, the M refers to the numerical order 
of the month beginning witb January as 1, 

MD,T, low-flow characteristics are the annual minimum D-day mean 

flow for T-year recurrence interval, in cubic feet per second, 
VD,T, flood volume characteristics are the annual maximum D-day 

mean flow for T-year recurrence interval, in cubic feet per 
second, and 

D50, discharge, in cubic feet per second, exceeded 50 percent of 
the time. 


Values for the 40 streamflow characteristics for each of the 20 
gaged basins used in. this investigation are available from the U.S. 
Geological Survey computer file, Streamflow/Basin Characteristics. 
Identical values were used in this investigation and in the Original 
streamflow data evaluation (Forrest and Walker, 1970). , 
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REGRESSION EQUATIONS 


Tables 3, 4, and 5 sunanarise the results of the multiple regression 
analyses. These analyses defined' mathematical equations of the form; 

log Y. = b(l) log X(l) + b(2) log. XC2)... + b(n) log X(n) + a + b(n+l) 
X(n+1) + b(n-b2) X(n+2)... + b(ra) X(ia) 
or its equivalent form; ' 

Y = X(l)^^^^X(2)^^^\..X(n)^^^\o l^Cn+l) X(n+1) + b(n+2) 

X(n+2)... + b(m) XCra)) 

where Y represents a streainflow characteristic, X(l) to X(m) represent 
'.>asin characteristics (where numbers in parentheses designate numeric 
order in a series), a represents the regression constant, b(l) to b(m) 
represent regression coefficients. In the program group, X(l) through 
X(n) were log transformed while X(n+1) through XCra) were not log trans- 
formed prior to calculating equations. Basin characteristics (independent 
variables) such as drainage area (A) and stream length (L) which have a 
larger range in values than other Independent variables for the 20 basins 
were log transformed. Those with a small range were not log transformed. 
No consideration was given in this study to the appropriateness of a 
regression model using solely logarithmic transforms, or to the reason- 
ableness of the exponents of characteristics included in derived equa- 
tions. Rather, a simple test to isolate effects of characteristics de- 
rived from Landsat by adding them to an existing data matrix and observing 
the changes in derived equations was performed. 
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In the tables the first column indicates the streamflow character- 
istic. In tables 3 and A, the second column presents the number of 
basins out of 20 that were used in the regression analysis. The next 
set of columns show the computed regression constant and the regression 
coefficients for those independent variables found to be statistically 
significant at the 95 percent level, and which significantly improved 
the accuracy of the equation. 
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Talile 3. — Control group equations relating streamElow character is tics to physiographic and cliiaatic 
characteristics of drainage basins as determined from maps and climatological records. 

,b(l) bC2) h(3)„hC4),rt(a + b(5)St + bC6)P + b(7)Si + b(8)P + b(9)I24,2 + b{10)Sn + b{ll)Tl + b(12)T7) 
Y — A 5 Ij “ 



Plow 

Number 

Regress- 







charac" 

of basins 

ion 


Regression coefficient 

for indicated basin characteristic 



teristlc 

used 

constant 







V 


a 

A 

.3 L E 

St 

F Si P 124,2 Sn 

T1 

T7 

P2 

20 

3.666 

0.806 

0.618 


-0. 963 



P5 

20 

10.065 

0.754 



-0.183 



PIO 

20 

10.931 

0.727 



-0.198 



P25 

IB 

"21.546 

0.447 

■* 


-0.007 


0.280 

i P50 

9 

-40.155 

0.499 





0.492 

u> 










18 

"1.581 

0.984 

0.188 


0.387 



SPA 

18 

-0.415 

1.019 






qio 

18 

-0.209 

0.908 






QU 

18 

-2. 003 

P. 961 



0,432 0.030 



Q12 

18 

-1.690 

1.022 

0.272 


0.364 



Q1 

18 

0-1. 074 

1.031 

0.131 


0.289 



Q2 

18 

-0,891 

1.014 


0.025 











Q3 

18 

-1.159 

1.018 





Q4 

18 

-0.763 

0.980 


C 

c 





18 

0.036 

0.973 


( 

1 

' trl 








1 

S Q 








1 

g 









St-* 









•^03 













Talile 3, — Coatrol group equations relating strearaflow characteristics to physiographic and cilmacic 

characteristics of drainage basins as determined from maps and climatological records-Continued. 
Y = + b(5)St + b(6)r + bC7)Sl + b(B)P + b(3)I24,2 + b(10)Sn + b(U)Tl + b(12)T7) 


T?loW 

cliarac- 

teristic 

. y 


Humber 
of basins 
used 


Regress- 

ion 

constant 
' a 


A 


Regression coefficient for Indicated basin characteristic 
L E St R Si R 124,2 


Sn 


T1 


T7 


I 

LJ 

o> 

1 


w 


o te 

i> § 

a® 

SB 

CO 

sa 


Q6 

18; 

-0.141 

0.987 

Q7 

18 

“0.162 

0.917 

QB 

18 

-0.078 

1.028 

Q9 

18 

-2.866 

0.860 

SDIO 

18 

-3-457 

0.951 

SDH 

18 

-0.318 

. 1.069 

SD12 

16 

-0.308 

: 1.116 

SDl 

18 

-0.185’ 

1,089 

SD2 

18 

0.228 

1.083 

SD3 

18 

-2.529 

1,049 

SD4 

18 

-0.162 

1.040 

SD5 

IS 

-0.261 

1.035 

SD6 

18 

-0.935 



SD7 

SD8 


18 

18 


-1.016 

-0.154 


0.973 


2,052 


0.598 0.043 


0.0B9 -0.466 


0.112 


-0.023 


0.042 


1.150 
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Table 3- — Control group equations relating streamElow characteristics to physiographic and climatic’ 

characteristics of drainage basins as determined from maps and climatological records-contlnued* 

b(l) bC2) h<3) b(4) (a + b(5)St + b(6)F + b (7)Si + b(8)P + b(9)I24»2 + b(10)Sn + b(ll)Tl + b(l2)T7) 
s L E 10 


Flow 

charac- 

teristic 

Y 

Humber 
of basins ■ 
used 

Regress- 

ion 

constant 

a 

A 

S 

•Regression coefficient for indicated basin characteristic 
L E St F Si P 124,2 

Sn 

T1 

T7 

SD9 

18 

-0.342 

0.914 


0.059 

-0.478 



M7,2 

.17 

. -1.389 



1.953 




117,10 

16 

-4.857 

• 


4.957 




H7,20 

15 

-24.151 



4.279 5.704 




V3,2 

10 

0.770 

1.080 






V3.25 

10 

1.001 

1.327 


-0.068 




V7,2 

17 

0.626 

1,071 






V7,1Q 

17 

0.919 

1,038 






V7,2S 

10 

0.900 

' 1.131 






E50 

11 ■ 

r.Q.451 



1.691 






y = A 


■ Table 4." — Experimental group equations relating streamflow characteristics to physiographic and climatic characteristics 
of drainage basins as determined from maps, climatological records, and Landsat imagery. 

b(l) b(2) b(3)^b(4)^^(a + b(S)5t + b(6^F + b(7)Si + b(8)P + b(9)I24,2 + b(10)Sn + b{U3Tl + b{12)T7 + b{13)Uan + bCl4)T3f + b(15)Urv + b(l6)Uw) 


Flow Number . Negress- 

charac- of basins ion 

teristic used constant- 

y a A- 


Regression coefficient for indicated basin characteristic 
S L. E St F Si P 124,2 Sn T1 T7 Uau tif Urv Ew 


P2 - 

20 

4.451 

0.717 

0.514 

-0.008 

-0.083 


P5 

20 

5.349 

0.964 0.570 



-0.092 


PIQ 

20 

-5,331 

0.569 

• 


-0.126 


P25 

18 

-21.546 

0.447 


-0.007 



P50 

9 

-40.155 

0.499 





QA 

la 

-1.581 

0.984 

0.1B8 



b.387 

SDA 

18 

-0.415 

1.019 





QIO 

18 

-0,209 

0.908 

■ 




Qll 

18 

0.049 

0.999 





Q12 

18 . 

-1.690 

1.022 

0.272 



0.364 

qi 

18 

“1.074 

1.031 

0.131 



0.289 

Q2 

18 

-0.891 

1. 014 



0,025 


Q3 

18 

-1.159 

1,018 



0.032 


Q4 

18 

-0.763 

0.980 


' t 


0.291 




0.009 -0.148 

-0.254 

0.153 -0.246 

0,280 

0.492 


-0.004 
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Table 4. — EiKperimeatal group equations relating streamflow characteristics to physiographic and climatic characteristics 

• V 

of drainage baains as determined from maps, climatological records, and landsat imagery -Contitiued. 


b(l) h(2) fa(3) bW) (a + b(5)St + b(6)F + bC7)S± + b{8)P + b(9)l24»2 + h<10)Sn + b(ll)Tl + b(12)T7 + b(l3)Uau + b(l4)Uf + b(lS)Drv + b(16)Uw) 
Y = A s L E 10 

Flaw 

charac- 

Futnber 
of basins 

Regress- 

ion 



' Regression coefficient for indicated basin characteristic 

ter Is tic 
Y 

used 

constant 

a 

A 

S 

L E St F Si P 124,2 Sn T1 17 TJau Uf Urv Uw 

Q5 

18 

0.036 

0.973 



Q6 

18 

-0.122 

1.048 


0.007 -0.011 

Q7 

18 

-0.076 

0.998 


-0.009 

Q8 

18 

-1.393 

1.024 


0.425 -0.003 

iO 

18 

-0.022 

0,964 


-0.008 0.154 

‘ SDIO 

18 

-3.457 

0.951 


0.112 

SDll 

18 

-0.318 

1.069 



SD12 

18 

-0.308 

1.116 . 



SDl 

18 

-0.140 

1.074 


0.077 

SD2 

18 

-0.268 

1.073 


0.004 

SD3 

18 

-2.529 

1.049 


0.089 -0.466 

SD4 

18 

-0.162 

1,040 



SD5 

18 

-0.261 

1. 035 

• 


SD6 

18 

-0.935 



2.052 

SD7 

18 

-0.819 

0.938 


0.007 

SDB 

18 

-0.154 

1.150 







Table A.— Experimeiital group equations relating ^streamflow characteristics to physiographic and climatic characteristics . 
of drainage baslna as determined from maps, climatological records, and landsat Imagery-Continued. 

hCl) hf2) b(3) b(A) (a + b(5)St V b(6)P + b(7)Sl + b(8)P + bC9)I2A,2 + b(10)Sn -i- bCll)Tl + b(12)T7 + b(13>Uau + b(14)Uf + b(15)Urv + b(l6)llH) 
Y = A S L E 10 


Flow Number Regress- 

charac— of basins Ion . 
teristic used constant 

, Y .. a 


Regression coefficient for indicated basin, characteristic 


T1 T7 


Uau Uf • Urv 


SD9 

18 

-0,342 

0.914 

0.059 

-0.478 


- 

M7,2 

17 

-95.603 


2.121 

• 

-0. 015 


M7,10 

16 

4.760 


5.543 0,005 


-0,008 


M7,20 

15 

-24.151 


4.279 

5.704 



V3,2 

10 

0.B50 

1,085 




t0.183 

V3,25 

10 

1,001 

1.327 

-0. 068 




TJ,2 

17 

0.691 

1.047 




-0,092 

V7,10 

V7,25 

17 

10 

0.919 
1. 048 

1-058 

1.148 

-0.002 

ss 

sg 


“0.097 

D50 

11 

-0.451 


1.691 






4 . 


Table 3.“-CbEiipar-i.aon of equatioas as to change in standard error of estimate of the streamflow character- 
istic resulting from use of Landsat-derived basin characteristics. 

S. E. : Standard error, in percent. 

AS* E. ; Improvement in standard error, in percent, caused by addition of Landsat-derived 


characteristics . 



Flow Independent variables S.E,, AS.E., 

charac- included in log units in percent^ in percent 

teristic 


Y Control group Experimental group Control Exper. Control Exper. 



P2. 

A.; 

E; 124,2 

A; 

E; 

F; P; Uf; Uw 

0.159 

0.092 

37.4 

' 21.2 

43.3 

T’5 


F 

A; 

s; 

P; Uw 

0.197 

0.134 

46.5 

31,5 

32.2 

Flo 


F 

A; 

F; 

T7; Uw 

0.205 

0.149 

49.0 

34-9 

28.8 

P25 

. A; 

F; T7 

A; 

F; 

T7 

0.154 

0.154 

36.3 

36.3 


P50 

A; 

T7 

A; 

T7 


0.229 

0.229 

51.1 

51,1 


QA 

A; 

E; 124,2 

A; 

E; 

124,2 

0.054 

0.054 

12.8 

12.8 


SDA 

A 


A 



0.088 

0.088 

20,1 

20.1 


QIO 

A 


A 



0.149 

0.149 

28.8 

28.8 


Qll 

A; 

124,2; Sn 

A; 

Uf 


0.098 

0.111 

23.0 

25.8 

12.2^ 

Q12 

A; 

E; 124,2 

A; 

E; 

124,2 

0.045 

0.045 

10.2 

10.2 



. A; 

E; 124,2 

A; 

E; 

124,2. 

0.047 

0.047 

10.8 

10.8 


Q2 

A; 

F 

A; 

P 


0.054 

0.054 

12.8 

12.8 


Q3 

A; 

P 

A; 

F 


0.062 

0.062 

14.3 

14.3 ' 



Table 5. — Comparison of equations as to change in standard error of estimate of the streamflow character 
istic resulting from use of Landsat-derived basin characteristics -Continued . 

S. E. : Standard error, in percent. ' 

A s. E.: Improvement in standard error, in percent, caused by addition of Landsat-derived 

characteristics. 


Floxj Independent variables S.E., • AS.E., 

charac- included in log units in percent;^ in percent 

teristic 


Y 

Control group 

Exp er imental g roup 

Control Exper. 

Control 

Exper. 


Q4 

A; 124,2 

A; 

124,2 

0.078 

0.078 

17.1 

.17.1 



A 

A 


0.120 

0.120 

28.0 

28.0 


Q6 

A 

A; 

P; Uf 

0.153 

0.097 

36.2 

22.8 

37.0 

Q7 

A 

A; 

tJf 

0.238 

0.168 

57.2 

39.6 

30.8 

Q8 

A 

A; 

124,2; Uf 

0.134 

0.106 

31,5 

24.6 

21.9 

Q9 

A; 124,2; Sn 

A; 

Uf; Uw 

0,136 

0.102 

31.9 

23.8 

25.4 

SDIO 

A; T1 

A; 

T1 

0.159 

0,159 

37.2 

37.2 


SDll - 

A ^ 

A 


0.165 

0.165 

39.0 

39.0 


SD12 

A 

A 


0.139 

0.139 

32.5 

32.5 


SDl 

A 

A; 

Uw 

0.083 

0.073 

19.3 

17.0 

11.9 

SD2 

Aj Sn 

A; 

Uf 

0.102 

0.101 

23.8 

23.5 . 

1.3 



Table 5. — Comparison of equations as to change in standard error of estimate of the streamflow character- 
istic resulting from use of Landsat-derived basin characteristics -Continued , 

S. E. : Standard error, in percent. ■ ' 

AS. E.: Improvement in standard error, *' percent, caused by addition of Landsat-derived 

characteristics. 


Flow 

charac- 

teristic 

Y 

Independent variables 
included 

Control group Experimental group 

S ,E. , 

in log units 
Control Exper. 

in percent^ 
Control Exper. 

As.e., 

in percent 

SD3 

A; P; 124,2 

A; P; 124,2 

0,084 

0.084 

19.4 

' 19.4 


SD4 

A . 

A 

0.088 

0.088 

20.1 

,20.1.’ 


SD5 

A 

A 

0.106 

0.106 

24.6 

24.6 


SD6 

I' 

L 

0.158 

0.158 

37.0 

37,0 


SD7 

A; Sn 

A; Uau 

0.225 

0.213 

54.0 

,50.8 

5.9 

SD8 

A 

A 

0.160 

0.160 

37.5 

37.5 


SD9 

A; E; Sn 

A; E; Sn . • 

0.158 

0.158 

39.5 

39.5 


M7,2 

L 

L; Uf 

0.482 

0.388 

136 

102 

25.0 

M7,10 

L- . 

'l; F; Uf 

1.220 

0.825 

No meaningful 

equation 

M7,20 

Lj 124,2 

L; 124,2 

1.264 

1.264 

No meaningful 

equation 

V3,2 


A; Uw ' 

0.127 

0.051 

29.5 

12.0 

59.3 


Table 5. — Cotnparison of equations as to change in standard error of estimate of the streamflow character 
istic resulting from use of Landsat-derived basin characteristics -continued . 

S. E. : Standard error, in percent. • 

As. E. : Improvement in standard error, in percent, caused by addition of Landsat-derived 

characteristics. . , 


Flow Independent variables S*.E. , A S.E. 

charac— included in log units in percent^ in percent 

teristic ‘ ■ ' ' 

Y Control group Experimental group Control Exper. Control Exper- 


V3,25 

A; St 

A; St 

0.112 

0.112 

26.1 ' 

26:1 


V7,2 

A 

A; Uw 

0.088 

0.073 

20.1 

17.0' ' 

15.4 

V7,10 

A 

A 

■ 0.114 

0.114 

26.4 

26.4 


V7,25 

A 

A; F; Uw 

0.094 

0.048 

21,9 

10.8 

50.7 

D50 

L 

L 

0.181 

0.181 

42.8 

42.8 


^Standard error, 

in approximate equivalent 

percent, 

calculated 

from standard 

error in 



logarithmic units (Thomas and Benson, 1970, p. 31). 

^Standard error increased by 12.2 percent; equation not improved. 



Table 5 suimnarizes the differences betT'/een independent variables 
included in the control and in the experimental equations. The standard 
error of each equation is given in logarithmic units and also in approX“ 
imate equivalent percentage. The percentages are actually the arithmetic 
averages of the plus and minus percentages about the mean calculated from 
the standard error in log units. Thus, an average standard error of 24.4 
percent represents a deviation of 27.3 percent on the plus side and 21-5 
percent on the minus side of the mean of the streamflow characteristic 
(Benson, 1962) . The last column shows the change in standard error where 
one or more variables derived from Landsat were included in the experimen- 
tal group of .equations. The change Is expressed as a percentage of the 
old (control group) standard error in log units. ' 

The value of the change is given for all equations computed except 
the 7-day, 10-year low flow, and the 7-day, 20-year low flow. The stan- 
dard errors for these equations were not significantly less than the 
standard deviatipns of the values of the dependent variables. 

An improvement in the acctiracy of the control group equations is 
considered substantial when there is at least a 10-percent reduction in 
the standard error of estimate of the control group equations by includ- 
ing basin characteristics derived from Landsat imagery. Table 5 shows 
that the standard errors of 15 out of 40 equations x^ere changed. Fourteen 
equations xrare improved and one equation xaas not improved. Among the 
14 equations improved, 12 x^ere improved by at least 10 percent. 
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OCGURRJSKCE OF VARIABLES 

IN EQUATIONS 





Table 5 shox-rs that variables derived 

from Land sat were included 

18 

times in the experimental group equations. 

, The most often used var- 


iables were forest and 

x-7ater as indicated 

beloxj: 




Streamflox<7 


Number of 

times that 


charact eristic 

Total Number 


indicated variable 


group 

of equations 


occurred 

in equations 



Uau 

Uf 

Urv 

Uw 

High 

10 

0 

1 

0 

6 

Average 

14 

0 

5 

0 

1 

Low 

3 

0 

2 

0 

0 

Variability 

13 

1 

1 

0 

1 

All characteristics 

40 

1 

9 

0 

8 


Forest, Uf, was included in 5 out o£ 14 equations for streainflow char- 
acteristics describing average flox<7* The coefficients (table 4) for 
forest are negative and imply an inverse relation between forest cover 
and mean monthly streainflow for. summer end early fall months in the 
study area. Forest \>ras also included in 2 out of 3 equations for low 
flows. In both cases the coefficients are also negative, implying an 
inverse relation between forest cover and annual minimum. 7-day mean 
flow. . 
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Water, Uw, was included in 6 out of 10 equations for high flows. 

In all 6 occurrences, the coefficients for water are negative, implying 
an inverse relation between v?ater and flood peaks and volumes with less 
than 25-year recurrence. Water as defined and measured in this investi- 
gation is considered to he a measure of surface water storage. 

Combined agricultural and urban land use, Uau, was included in only 
1 out of 13 equations for flow variability. Riparian vegetation, Urv, 
as defined and measured in this investigation was not included in any 
of the equations for strearaflow characteristics. 

Inferred hydrologic slgnific .nee was used in the initial selection 
of all basin characteristics. However, the basis for including any 
characteristic in an equation is primarily statistical. The inter-rela- 
tions between the basin characteristics along with the inability of 
the characteristics to completely describe a drainage basin makes tenuous 
any assertions about the physical effects of the basin characteristics 
on streamflow. Despite the inability of the relations to describe the 
fundamental causes of streamflow variations, the basin characteristics 
frequently included in the aquations are numerical measures that are 
related to the flow variations. 



Regression analyses were performed with hoth measures of forest, 

P (map“derived) and Uf (Landsat-derived) , and surface water storage, St 
(map-derived) and Uw (Landsat-derived), in the data matrix. This exped- 
ient was expected to result in one measure being replaced by the other 
where significance was obtained. In 3 cases (P2, Q6, and ^^/,10) how- 
ever, the inclusion of Uf generated the inclusion of P. In all 3 cases 
the exponents of P and Uf are of opposite sign. Because the simple corre- 
lation coefficient of F versus Uf is 0.82, inclusion of hoth in the 
derived equation tends to cancel their effect and makes the computed stan- 
dard error suspect. A conclusion which is still valid for the 3 cases, 
however, is that Uf is a more powerful characteristic than P; other^rlse, 

P xrould have appeared in the control group equations, ^-Jhere only Uf 
appears, the standard error improvement criterion is valid as a demon- 

stratable effect of using Landsat-derived data. 

Concurrent appearance of St and TJvj did not occur. The simple" corre- 
lation coefficient of St versus Uw is 0.02, and the two deteimtinations 
were completely independent, 

P and St were included in only one control group equation each* 

In contrast to this, Uf and Uw were included in 9 and 8 experimental 
group equations, respectively. Accordingly, Uf may be a more reliable 
measure of forest cover in the study area than F, and Uw may he a more 
reliable measure of basin storage than St, If so, the frequent inclu- 
sion of Uf and Uw in the experimental group equations is consistent 
with the frequent inclusion of F and St in equations for the Potomac 
River basin immediately west of the study area (Thomas and Benson, 

1970). 
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In the Potomac River basin study, however, F was included in 11 out 
of 14 high flow equations and in only 2 out of 15 average flox? equations. 
These results appear inconsistent with results of this investigation un- 
less differences in the average physiography of the two study areas are 
considered. In the flat, sandy tcrrc.n of the Delmarva Peninsula, 
forest should be more closely related to evapotranspiration and therefore 
to low- flox^s as X'/ell as mean monthly flox^s for summer and fall months. 

In the hilly, rocky terrain of the Potomac River basin, forest should be 
more closely related to steep untillable slopes and therefore to flood 

f lOXxTS . 

The basin characteristic riparian vegetation, Urv, was not accurate- 
ly isolated during image analysis, and it contains considerable upland 
vegetation. This may be the principal reason why riparian vegetation 
was not included in any of the experimental group equations. 
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SUMMARY MD CONCLUSIONS 


This study tested the usefulness of basin characteristics derived 
from Landsat imagery for improving equations used to estimate streamflow 
characteristics i The Delmarva Peninsula of Maryland, Delaware, and 
Virginia is a study area representative of rural land use in areas of 
low topographic relief on the humid east coast of the United States. 

Basin characteristics derived from Landsat imagery, especially forest (Uf) 
and x^ater CUtj), are representative of characteristics x^hich most readily 
can he extracted from Landsat imagery by multispectral image analysis sys- 
tems employing film density-discrimination techniques. Other hydrolog- 
ically-significant characteristics may be extracted as remote sensing 
technology improves. 

The basin characteristics derived from maps and climatological records 
and used in the control group of equations computed by the multiple- 
regression program xcere the same characteristics used in the original 
streamfloxT data program evaluation (Forrest and Walker, 1970) . They cover 
a x^ide range of characteristics found to be significant in several pre- 
vious studies and include, contributing drainage area, main channel slope, 
stream length, mean basin elevation, area of lakes, ponds, and sx^amps, 
forest area, soil infiltration index, mean annual precipitation, precipita- 
tion intensity, snox'/fall, and mean minimum January and mean maximum July 
temperature. 
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The streamflow characteristics i,^j.ived from the records of daily 
discharge of 20 gaged basins in the study area are representative of 
the full range in flow conditions and include all of those commonly 
used for design or planning purposes. They include annual flood peaks 
x^ith recurrence intervals of 2, 5, 10, 25, and 50 years, mean annual 
discharge, standard deviation of the mean annual discharge, mean monthly 
discharges, standard deviation of the mean monthly discharges, low-flow 
characteristics, flood volume characteristics, and the discharge equalled 
or exceeded 50 percent of the time. 

These streamflow characteristics were related to the basin char- 
acteristics of the corresponding 20 drainage basins by a technique of 
multiple regression using a digital computer. A control group of equa- 
tions was computed using basin characteristics derived from maps and 
climatological records. An experimental group of equations was computed 
using basin characteristics derived from Landsat imagery as well as from 
maps and climatological records. The standard error of estimate of the 
two groups of equations was compared to see if any reduction in standard 
error could be considered a substantial improvement upon the original 
equations. 

Based on a reduction in standard error of estimate equal to or 
greater than 10 percent, the equations for 12 streamflow characteristics 
were substantially improved by adding to the analyses basin characteristics 
derived from Landsat imagery. These improvements are summarized in 
table 6. 

■ _5i„- ■■ 



Table 6. — Twelve streamflow characteristics and corresponding standard 
errors reduced by at least 10 percent. 

Standard error, 




in 

percent 

Change in 




Experi- 

standard error 



Control 

mental 

in percent of 


Streamflow characteristic 

group 

group 

control group 

P2, 

annual flood peak of 2-year 





recurrence 

37.4 

21.2 

43.3 

P5, 

annual flood peak of 5™year 





recurrence 

46.5 

31.5 

32.2 

PIO 

, annu^'-l flood peak of 10-year 





recurrence 

49.0 

34.9 

28.8 

Q6, 

mean monthly flow for June 

36.2 

22.8 

37.0 

Q7, 

tt It Ti .1 

57.2 

39.6 

30.8 

QS, 

" " ” August 

31.5 

24.6 

21.9 

Q9, 

" . " " September 

31.9 

23.8 

25'. 4 

SDl 

, standard deviation of mean 





flow for January 

19.3 

17.0 

11.9 

M7, 

2, annual minimum 7~day mean 





floi'j for 2-year recurrence 

136 

102 

25.0 

V3, 

2, annual maximum 3-day mean 





flow for 2-year recurrence 

29.5 

12.0 

59.3 

V7, 

2, annual maximum 7-day mean 





flow for 2-year recurrence 

20.1 

17.0 

15.4 


V7,25, annual maximuni 7-day mean 

flow for 25-year recurrence 21.9 10.8 50.7 
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Improvements occurred In all flow regimes. The hasin characteristics 
derived from Landsat imagery that were included most frequently in the 
experimental group equations are forest, Uf, and water, Uw. It is possible 
that these tw characteristics are more reliable measures of forest cover 
and basin storage than F and St used in the previous streamflow data eval- 
uation for Maryland and Delaware (Forrest and Walker, 1970). 

By comparing the reduced standard error of each streamflow character- 
istic with the accuracy goal for tliat characteristic, it is possible to 
judge the degree to which some goals can be achieved by using basin char- 
acteristics derivei. from landsat imagery. In those cases where the goals 
can thus be achieved, the streamflow data collection effort can be re- 
directed to areas of higher priority. 
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